Third and Fourth Order Phase Transitions: 
Exact Solution for the Ising Model on the Cayley Tree 
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An exact analytical derivation is presented, showing that the Ising model on the Cayley tree 
exhibits a line of third order phase transition points, between temperatures T2 = 2/c^ 1 Jln(v2 + 1) 
and Tbp = kg 1 Jln(3), and a line of fourth order phase transitions between Tbp and 00, where ks 
is the Boltzmann constant, and J is the nearest-neighbor interaction parameter. 
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Ever since the pioneering works on phase transitions 
and critical phenomena, it has been clear that phase tran- 
sitions of higher order are conceptually possible, but, to 
the best of our knowledge, up to date there has been no 
rigorous proof of existence of a single system exhibiting 
a phase transition of finite order higher then two. 

One notable attempt in this direction is the work of 
Muller-Hartmann and Zittartz interpreting the se- 
ries expansion of the free energy of the Ising model on 
the Cayley tree in the limit of small field H — ► + , 
in terms of phase transitions of all even orders (from 
two to infinity). More precisely, their claim was that 
the susceptibility of order 2£ (i.e. the 2£-th deriva- 
tive of the free energy with respect to field), diverges 
at temperature T 2 e = fc^Jlnfl + 2/(B 1 ~ 1 / 2 < ! - 1)] for 
£ = 1,2,. ..,00, where B is the tree branching num- 
ber (coordination number minus 1). Thus, the zero 
field susceptibility was found to diverge from T — to 
T-2 = kg 1 J\n[(yB + \)/(\/B~ — 1)], representing a line 
of second order phase transitions (with respect to field), 
fourth order susceptibility was found to diverge between 
T 2 andT 4 = k B x Jln^B 3 / 4 + 1)/(S 3 / 4 - 1)], representing 
a line of fourth order phase transitions, while the infinite 
order phase transition was found to occur at the Bcthe- 
Pierls temperature T B p = k^Jln [(B + 1)/(B - 1)]. 

It is shown in this work through an exact analytical cal- 
culation of the recursion relations for the derivatives of 
the free energy, that the exact expressions contain terms 
previously unnoticed in the approximate analysis. In- 
terpratation of these results requires some special con- 
sideration of symmetry breaking, the classical approach 
not being adequate for the system at hand. It is found 
that the third order derivative diverges between T 2 and 
Tbp (all the lower derivatives being finite), representing 



a line of third order phase transitions, while the fourth 
order derivative diverges between Tp p and 00 (where all 
lower derivatives are finite), representing a fourth order 
phase transition line. 

We consider the nearest-neighbor Ising model with the 
Hamiltonian 
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-J SiSj — H^^Sj 
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where J is the coupling constant, H is the external mag- 
netic field, Si = ±1 is the spin at site i, and (nn) de- 
notes summation over the nearest- neighbor pairs. For 
simplicity, hereafter we consider only the tree with B = 2, 
while the analysis which follows may be directly gener- 
alized for arbitrary tree branching number. Following 
Eggarter 0], we further consider systems situated on a 
single n-generation branch of a Cayley tree, composed 
of two [n — l)-generation branches connected to a sin- 
gle initial site. Thus, the n-generation branch consists 
of N n = 2" +1 — 1 spins, the 0-generation branch being 
a single spin. The exact recursion relations for the par- 
tial partition functions of any two consecutive generation 
branches are easily derived [2| to be 
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where and Z~ denote the partition functions re- 
stricted by fixing the initial spin (connecting the two 
n-generation branches) into the {+} and {— } position, 
respectively, and where we have used the notation x = 
exp(f3J) and y = exp((3H), with (3 = 1/kpT denoting 
the reciprocal of the product of the Boltzmann constant 
ks and the temperature T. The usual approach of at- 
tempting to establish the field dependent expression for 
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the partition function in the thermodynamic limit, then 
finding its field derivatives, and finally taking the limit 
H — ► yields only approximate solutions for the zero- 
field susceptibility U SO] • In an earlier work |j| , current 
authors have derived an exact expression for the zero-field 
magnetization and susceptibility, by using the strategy of 
finding the recursion relations for the field derivatives of 
the partition function, taking the limit H — ► 0, and 
only then performing the actual iterations to reach the 
thermodynamic limit. Here we extend this approach to 
find higher derivatives of the partition function. 

Equation J5J) can be formally differentiated with re- 
spect to field, to find the recursion relations for the field 
derivatives of the partition function. Up to the forth 
derivative we find the following recursion relations 



Starting from a single spin (0-th generation branch), 
for which we have d k Z^ /dh k = (±l) fe y ±1 , it is straight- 
forward to show by mathematical induction, using (12l(il) . 
that for zero field [y = 1) the symmetry equations 
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hold for a branch of arbitrary generation n. It then 
follows from I12I7II that the moments S n = -jr^§j^, 



r = t , 

n ~ Z+ 9h 2 ' — z+ dh.3 ' 

the recursion relations 



and V n 
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satisfy 



92 n+l = ±1 / . _p, ! - 

dh u '~ 



d 2 z 



n+l 



dh 2 



^ Z n+1 

dh 3 



0±T^1± 



2 pu± r 



= y 



±i 



r r o±N 2 ±4r o±pi± 

I n ) n *• n 



2 r. 



i±-\ 



(3) 



(4) 



,±i 



±(r°±) 2 + 6r°±r^±6(r^) 2 ± 



± gp0±p2± , gpl±p2± 



ar^r^l, (5) 



Sn+l — 1 + 2 t S n , 



%i+l — 1 
Un+l = 1 



V 



n+l 



■4tS n + 2t 2 S 2 +2T n , 
6tS n + 6t 2 S 2 + 6T„+ 
+6 t S n % + 2 1 U n 
8 t S n + 12 t 2 S 2 + 12T„ 
+24 tS n T n + 8t U n 
6 7^ + 8 t 2 5„ Z4 + 2 V n 



1 



So = 1, 
T = 1, 
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where we have used notation h = [3H and 
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where t = tanh(/3J). Relations ijHJ can be iterated (by 
summing the geometric series) to yield closed-form ex- 
pressions for S n , T n ,U n and V n for arbitrary tree gener- 
ation n, and we find 
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Before proceeding with the analysis of the above ex- 
act expressions, a word is due on symmetry breaking. 
It has been commonly accepted for this model that all 
odd derivatives of the free energy with respect to field 
are identically zero in zero field because of symmetry, 
and only even derivatives have been analyzed [ij, la, |7j (in 
fact, equations similar to (|10|l and (|12|l for the second and 
fourth derivatives have been derived in references 0, Q 
by neglecting odd correlations). However, this is true for 
any lattice in strictly zero field, as every spin configu- 
ration has a mirror image (obtained by flipping all the 
spins) with exactly the same energy, and inverted sign 
of the configurational magnetization ^ SV The usual 
procedure of breaking the symmetry by retaining an in- 
finitezimal positive field while taking the thermodynamic 
limit, and only then taking the zero field limit, is not suit- 
able in the present case because the derivatives of the free 
energy diverge in wide temperature regions (rather than 
just in a single critical point). Here it seems more appro- 
priate to break the symmetry by applying an infinitely 
strong local field of infinitezimal range, which can be im- 
plemented by restricting a single spin into one of the two 
possible orientations, while considering all the possible 
orientations of all the other spins. In a recent work |l0j . 
the present authors have analyzed the effect of restrict- 
ing a single spin in the {+} orientation on magnetization 
(first derivative of the free energy with respect to field), 
for the current system. It was found that symmetry is 
indeed broken by fixing an arbitrary (surface, bulk or cen- 
ter) spin. While without this restriction magnetization is 
identically zero in strictly zero field for arbitrary system 
size, it was shown that fixing any spin leads to magnetic 
ordering of extremely large systems, in a wide tempera- 
ture range (even if magnetization does go to zero in the 
thermodynamic limit for all nonzero temperatures). 

In the rest of this paper, we shall therefore adopt the 



strategy of breaking the symmetry by fixing a single (cen- 
tral) spin in the {+} (upward) orientation while taking 
the thermodynamic limit, and we shall henceforth use 
the term ^-th order susceptibility for the quantities 
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for finite size systems, and 
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for the thermodynamic limit. 

The restricted magnetization, < m > + = is found 
5, 8] to be zero for all nonzero temperatures in the 
thermodynamic limit, (even if it retains nonzero val- 
ues in a wide temperature range for systems far ex- 
ceeding in number of particles the observable Universe 
[loj|'). while 2-nd order susceptibility Xn^ diverges below 
T 2 = 2fc B 1 Jln(\/2 + 1). It was also recently shown [9j 
that in the thermodynamic limit the divergence of x^ 2 ^ 
in the vicinity and at Ti is extremely weak, with critical 

(2) 

exponent 7 = 0, where susceptibility of a finite tree \n 
diverges proportionally to the three generation level n, as 
n — > 00 . In the rest of this paper, we analyze the third 
and fourth order susceptibility given by Xn" 1 = f3 2 U n /N n 
and Xn^ = f3 3 V n /N n , respectively. 

Equations I10ll2f> are deceptive in the sense that at 
first glance they suggest divergence of second, third and 
fourth order susceptibilities for arbitrary tree generation 
level n, at points 2t 2 = 1 and 2t = 1 (corresponding to 
T2 and Tbp, respectively), while in addition, the fourth 
order susceptibility seems to diverge (irrespective of n) at 
4< 3 = 1 and 8i 4 = 1 (the last corresponding to the Muller- 
Hartman and Zittartz temperature T4 [l|). In fact, the 
exact expression corresponding to i|12|) . but obtained by 
neglecting odd correlations 6], indeed does diverge at 
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TABLE I: Leading terms (for large n) of the series expansion 
of susceptibilities xS? , Xn' and Xn' i arround points 2t = 1 
and 2t 2 — 1, demonstrating that expressions H10I12H do not 
diverge at these points for any finite n (see text for details). 
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TABLE II: Leading terms of Xn 5 arround 4t 3 = 1 and 8t 4 = 1. 
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FIG. 1: Scaled second, third and fourth-order zero-field sus- 
ceptibilities, calculated using formulas (1101121 . for several sys- 
tem sizes n = 64, 128, 256 and 1024. The dotted vertical lines 
indicate the points T 2 and Tb p , where the transition changes 
from second to third, and from third to fourth order, respec- 
tively. 



8t 4 = 1. On the other hand, a more detailed analysis 
of (|10I12(I shows that the interplay between the individ- 
ual terms cancels out these aparent singularities at the 
mentioned temperature points. In particular, upon se- 
ries expansion of the exact expressions (|10I12|) arround 
these points (from either side), it is found that there is 
no divergence for any finite n. We give the leading terms 
of these expansions for large n in Tabs. HI and ITU 

Consequently, divergence of higher order susceptibili- 
ties in the thermodynamic limit in different temperature 
regions is caused only by diverging generation level n, 
and the finite size scaling is accomplished simply by the 
formula lnx„/n. In Fig. 1 we show the scaled curves of 
higher order derivatives, obtained using formulas l|lUll2fl . 
for several system sizes. 

It is seen from Fig. 1 that there are three distinct tem- 
perature regions, where each temperature represents a 
point of second, third, or fourth order phase transition. 
Between zero and T 2 the second derivative diverges (the 
first derivative being finite), and we have a second or- 
der phase transition line. Between T2 and Tbp the third 
derivative diverges (first and second derivatives being fi- 
nite), representing a line of third order phase transitions, 
while the fourth order derivative diverges between Tbp 
and 00 (where all lower derivatives are finite), represent- 
ing a fourth order phase transition line. The explicit 
expressions for the limiting curves, corresponding to the 
finite size system curves shown in Fig. 1, are obtained by 

taking the limit k 1 = lim„ ^ ln^^/n , and the results 

obtained using formulas (|10I12|) for different temperature 
regions are sistematized in Tab. IIIII 

In conclusion, after decades of continuous interest, the 
Ising model on the Cayley tree continues to furnish new 
insights into critical phenomena. Being a highly unphys- 
ical system (with its infinite dimension and finite order of 
ramification), Cayley tree may turn out a unique struc- 



TABLE III: Scaled second, third and fourth order suscepti- 
bility, in different temperature regions. 

0<T<T 2 T 2 <T<T B p T>T B p 

(l>t>l/\/2) (l/v / 2>t>l/2) (*<l/2) 

1P> ln(2t 2 ) 

k (3) ln(4t 3 ) ln(2t) 

k (4) ln(8t 4 ) In 2 In 2 



ture where phase transitions of order higher then two 
indeed do exist, and the latter may prove to be of only 
academic interest. On the other hand, it is possible that 
current findings may turn out relevant for interpretation 
of experimental data on finite size branching structures, 
of the real observable physical world. 
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